PT-Symmetric Matrix Quantum Mechanics 
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Recently developed methods for PT-symmetric models are applied to quantum-mechanical matrix 
models. We consider in detail the case of potentials of the form V = -(g/N p/2 ~ 1 )Tr (iM) p and show 
how the calculation of all singlet wave functions can be reduced to solving a one-dimensional PT- 
symmetric model. The large- iV limit of this class of models exists, and properties of the lowest-lying 
singlet state can be computed using WKB. For p = 3,4, the energy of this state for small values of 
N appears to show rapid convergence to the large- N limit. For the special case of p = 4, we extend 
recent work on the — gx 4 potential to the matrix model: we show that the PT-symmetric matrix 
model is equivalent to a hermitian matrix model with a potential proportional to + (4g/N)TrTl 4 . 
However, this hermitian equivalent model includes an anomaly term h^2g / NTrli. In the large- N 
limit, the anomaly term does not contribute at leading order to the properties of singlet states. 



I. INTRODUCTION 

Matrix models appear in many contexts in modern theoretical physics, with applications ranging from condensed 
matter physics to string theory. Interest in the large- N limit of matrix models was strongly motivated by work on 
the large-iV c limit of QCD but interest today is much wider. For example, Hermitian matrix quantum mechanics 
leads to a construction of two-dimensional quantum gravity coupled to c = 1 matter 0] ■ 

We will show below that the matrix techniques pioneered in [3] for Hermitian matrix quantum mechanics can be 
extended to PT-symmetric matrix quantum mechanics, where the matrices are normal but not necessarily Hermitian. 
The large- N limit can then be taken in PT-symmetric matrix theories just as in the Hermitian case. Quantities of 
interest such as the scaled ground state energy and scaled moments can be calculated using WKB methods. In the 
special case of a quartic potential with the "wrong" sign, we prove using functional integration for all values of N 
that the PT-symmetric model is equivalent to a hermitian matrix model with an anomaly, as in the one-component 
case [HQ. Interestingly, the anomaly vanishes to leading order in the large- N limit. 
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II. FORMALISM 

The solution for all N of the quantum mechanics problem associated with the Euclidean Lagrangian 

2 \ dt ) N 



(1) 



where M is an N x N Hermitian matrix was first given by Brezin et al. [3j . The ground state ip is a symmetric function 
of the eigenvalues Xj of M. The antisymmetric wave function 4> defined by 
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satisfies the Schrodinger equation 
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(j) = N 2 E (0) (j) 



(2) 



(3) 



where E^°'is the ground state energy scaled for the large-./V limit. This equation separates into N individual 
Schrodinger equations, one for each eigenvalue, and the antisymmetry of cf> determines N 2 E^ as the sum of the 
N lowest eigenvalues. 
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Here we solve the corresponding problem where the potential term is PT-symmetric but not Hermitian. As shown 
by Bender and Boettcher [6j], the one- variable problem may be solved by extending the coordinate variable into the 
complex plane. This implies that for PT-symmetric matrix problems, we must analytically continue the eigenvalues 
of M into the complex plane, and in general M will be normal rather than Hermitian. We consider the Euclidean 
Lagrangian 

L= i Tr (^i) 2 _9^ Tr(iM)P (4) 
2 V dt J Np/ 2 - 1 K ' w 

with g > 0. Making the substitution M — > UAU + , with U unitary and A diagonal, we can write L as 



where 




In the analysis of conventional matrix models by Brezin et al., a variational argument shows that the ground state 
is a singlet, with no dependence on U. Because the Aj's are in general complex for PT-symmetric theories, this 
argument does not apply. However, in two cases we can prove that the ground state is indeed a singlet: for p = 2, 
which is trivial, and for p = 4, where the explicit equivalence with a hermitian matrix model proven below can be 
used. Henceforth, we will assume that the ground state is a singlet, but our results will apply in any case to the 
lowest-energy singlet state. 

We have now reduced the problem of finding the ground state to the problem of solving for the first N states of 
the single-variable Hamiltonian 

This Hamiltonian is PT-symmetric but in general not Hermitian. The case p — 2 is the simple harmonic oscillator. 
For p > 2, the Schrodinger equation associated with each eigenvalue may be continued into the complex plane as 
explained in [g|. We exclude the case p < 2, where PT symmetry is spontaneously broken and the eigenvalues of H 
are no longer real. 



III. GROUND STATE PROPERTIES 



As with Hermitian matrix models, the ground state energy is the sum of the first N eigenenergies of the Hamiltonian 
H. In the large N limit, this sum may be calculated using WKB. A novelty of WKB for PT-symmetric models is the 
extension of classical paths into the complex plane. This topic has been treated extensively in @, 0|- 

We define the Fermi energy Ep as the energy of the iV'th state 



N 



^ J dpd\9[E F -H(p,X)} (8) 



where the path of integration must be a closed, classical path in the complex p — A plane. In order to construct the 
large- N limit, we perform the rescaling p — > \ Np and A — > y/NX yielding 

H sc (p,X) = ±p 2 -g(tXr (9) 

where the scaled Hamiltonian H sc is related to H by H = NH SC . We introduce a rescaled Fermi energy cp given by 
Ep = Nep, which is implicitly defined by 

1 = ^ f dpdX9[e F -H sc (p,X)]. (10) 
After carrying out the integration over p, we have 

1 = ~J dX^2e F + 2g {iXfO [e F + g {iXfj (11) 
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where the contour of integration is taken along a path between the turning points which are the analytic continuation 
of the turning points at p = 2. This equation determines ep as a function of g. 
We define a scaled ground state energy E^°> by 

N-l 
fc=0 

The WKB result for the sum of the energies less than Ep can be written as 

N ~ l n 2 r 

J2 Ek = Y~ / dpd\H sc {p, X)6 [e F - H sc {p, A)] (13) 
fc=o n J 

so that in the large- N limit is given by 

E ™ = hj dpdX Hscip ' x)d [eF ~ Hsc{jpi A)] (14) 

The integration over p is facilitated by using equation (|10[) to insert a factor of giving 

£?<°) = eF _ J- / dpdA [ eF - ff sc ( p , A)] [e F - ff sc (p, A)] . (15) 



2tt . 

The integral over p then yields 

B W = e F - i- / dA [2e F + 2 3 (zA) p ] 3/2 [«* + (i\f] . (16) 



The turning points in the complex A plane are 



X- = [j) e^Wa-Vrt (17) 



A + = f^) 1/P e — 1 



We integrate A along a two-segment, straight-line path connecting the two turning points via the origin 
equation (JTDJ) for ep, we find 



tt\p/ r(3/2 + l/p) 



(!) 



K sin (tt/p) T(l + 1/p) 
and solving (Tl6|) for the scaled ground state energy we have 
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(18) 
Solving 

(19) 



3p + 2 3p+2 



tt^p/ r(3/2 + i/p) 

2, 



sin (?r/p) T(l + 1/p) 



2p 



(20) 



For p = 2, this evaluates to E^ — y/g/2 , in agreement with the explicit result for the harmonic oscillator. 

It is very interesting to compare the large- N result with results for finite N. The low-lying eigenvalues for the 
Hamiltonian p 2 — (ix) p have been calculated by Bender and Boettcher in Q for the cases p = 3 and p = 4; the case 
p = 2 is trivial. We can use their results by noting that the eigenvalues of our Hamiltonian H are related to theirs by 



Ej = 



fJ 



2/(p+2) 



Tp- 

2p/(p+2)jy(p~2)/(p+2) 3 



BB 



(21) 



Results for p = 3 and 4 and small values of N are compared with the large- N limit in Table 1. The energies for finite 
values of N rapidly approach the N — > oo limit. The approach to the limit appears monotonic in both cases, but with 
opposite sign. 
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N 


P=3 


p=4 


1 


0.762852 


0.930546 


2 


0.756058 


0.935067 


3 


0.75486 


0.935846 


4 


0.754443 


0.936115 


5 


0.754251 


0.936239 


6 


0.754147 


0.936306 


7 


0.754084 


0.936347 


8 


0.754043 


0.936372 


oo 


0.753991 


0.936458 



TABLE I: The scaled ground state energy i% at g = 1 for p = 3 and p — 4. 



The expected value of (TrM) for large N is given by 



N-l 

(TrM) = < A i> = o- / d P dA A0 - # (P. A )l • ( 22 ) 
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Calculations of higher moments (TrM n ) are carried out in the same manner. Upon rescaling, we find that (TrM) 
grows as TV 3 / 2 , and the scaled expectation value is given by 
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(TrM) 



±- J dpdXXe [e F -H sc (p,X)\ 



which reduces to 



H = i J d\ X^2e F + 2g (iX) p 8 [2e F + 2g (iXf] 
Using the same two-segment straight line path as before, we find that 
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cos l — 



r(3/2 + 1/p) 



r(i + i/p) 



p+2 r(i + 2/p) 

r(3/2 + 2/p)' 



(23) 



(24) 



(25) 



For p = 2, /i = 0, as expected for a harmonic oscillator. For p > 2, the expectation value /i is imaginary because (Xj) 
for each eigenstate of the reduced problem is imaginary Q. For p — 3, fi — — 0.52006i. For p = 4, = —0.772539?. In 
the limit p — + oo, /i goes to -i. This behavior is easy to understand, because in this limit, the turning points become 
degenerate at — i. 



IV. SPECIAL CASE OF TrM 4 



For the case of a TrM 4 interaction, we can explicitly exhibit the equivalence of the PT-symmetric matrix model 
with a conventional Hermitian quantum mechanical system. As in the single-variable case, there is a parity-violating 
anomaly, in the form of an extra term in the Hermitian form of the Hamiltonian, proportional to H. We show below 
that the anomaly term does not contribute at leading order in the large- N limit. 

The derivation of the equivalence closely follows the path integral derivation for the single- variable case 0, [1[ . The 
Euclidean Lagrangian is 



L = \ Tr (T> +^TrM*-fTrUi 



and the path integral expression for the partition function is 



Z = J [dM]exp(- J dtlA. (27) 



Motivated by the case of a single variable, we make the substitution 

M = -2Wl + iH 



(28) 



where H is an Hermitian matrix. Because M and H are simultaneously diagonalizable, this transformation is tanta- 
mount to the relation 



\j = — 2i*Jl + ihj 

between the eigenvalues of M and the eigenvalues hj of H. The change of variables induces a measure factor 



[dM] 



[dH] 



Det[VTTiH] 

where the functional determinant depends only on the eigenvalues of H. The Lagrangian becomes 

= l T {dH/dtf _ 2 _ _g 2 

2 1 + TV 



(29) 



(30) 



(31) 



at the classical level. However, following [fj, we note that in the matrix case the change of variables introduces an 
extra term in the potential of the form 



which can be written as 



The partition function is now 
[dH] 

-exp 



AV : 



-E— 

32 ^ 1 + 
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(dh 3 \ 
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_dhj 


\d\J_ 



1 2 



ill, 



= Tr 

32 



1 + iH 



dt 



det [VTTiH] 

We introduce a hermitian matrix-valued field II using the identity 



1 (dH dt) 2 , m , , j m , N9 1 / 1 
-Tr^ - 2m 2 Tr(l + iH) -Tr 1 + iH) 2 Tr 

2 1 + iH v ' N y ' 32 \l + iH 



1 



■/ m "f/ 4IV 



2 v ; \ 1 + iH 



det [VI + 

Dropping and adding appropriate total derivatives and integrating by parts yields 

Z = j [dH]{dn}exp{- j dtTr -2m 2 (l + iH) - 16-^(1 + iH) 2 + i (1 + iH) n 2 + 11(1 + iH) 

The integration over H is Gaussian, and the shift H — > + i gives 

JV 



-n 



Z = / [dnjexp^ - / eftTr 



, ri 2 - 2m 2 n 2 + in 4 ) - in 

64g \ 4/4 



After the rescaling II — > y/32g/NH we have finally 

Z = y [dTI] exp j - y dtTr i (n 2 - 2m 2 n 2 ) + ^n 4 - y/2g/NU 
We have now proven the equivalence of the PT-symmetric matrix model defined by 



1 fdMY 1 
L = -Tr 1 4- - 

2 V dt 



m 2 TrM 2 - —TrM 4 
N 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 
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to the conventional quantum mechanics matrix model given by 




Trn 4 . 



(40) 



This equivalence implies that the energy eigenvalues of the corresponding Hamiltonians are the same. This could also 
be proven using the single-variable equivalence for the special case of singlet states, but the functional integral proof 
encompasses both singlet and non-singlet states at once. The equivalence of these two models also allows for an easy 
proof of the singlet nature of the ground state. Standard variational arguments show that the ground state of the 
Hcrmitian form is a singlet. The direct quantum mechanical equivalence of the single- variable case is then sufficient 
to prove that the ground state of the PT-symmetric form is also a singlet. 

As in the single- variable case, there is a linear term of order h appearing in the Lagrangian and Hamiltonian of the 
Hermitian form of the model. This term represents a quantum mechanical anomaly special to the TrM model. To 
determine the fate of the anomaly in the large- TV limit, we construct the scaled Hamiltonian of the Hermitian form 
in exactly the same way as for the PT-symmetric form. It is given by 



indicating that the effect of the anomaly is absent in leading order of the large- N expansion. One easily checks for 
the m = case that the Hermitian form without the linear term reproduces the PT-symmetric prediction for at 
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